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Abstract 

We study a model consisting of ./V nonlinear oscillators with global periodic 
coupling and local multiplicative and additive noises. The model was shown to 
undergo a nonequilibrium phase transition towards a broken-symmetry phase 
exhibiting noise-induced "ratchet" behavior. A previous study [10] focused 
on the relationship between the character of the hysteresis loop, the number 
of "homogeneous" mean-field solutions and the shape of the stationary mean- 
field probability distribution function. Here we show -as suggested by the 
absence of stable solutions when the load force is beyond a critical value- 
the existence of a limit cycle induced by both: multiplicative noise and global 
periodic coupling. 
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I. INTRODUCTION 



The study of dynamical systems has shown that limit cycles are ubiquitous in a wide 
range of physical applications [1,2]. From a physicist's point of view, limit cycles are thought 
of as a way to balance the in- and out- energy flows. Even when those flows are not 
oscillatory in time, a system's oscillatory motion can occur equalizing such flows over one 
period. A nice pedagogical example of such a process, based on a perturbative analysis of 
the nonlinear van der Pol oscillator, can be found in [3]. As is well known, limit cycles 
are robust -structurally stable under small perturbations- attractors in dissipative systems 
without external oscillations [1,2]. Usually, limit cycles arise in dynamical systems described 
by ordinary differential equations (ODE) [1,2], but there are several examples where such 
kind of cycles also arise in partial differential equations (PDE) or "extended systems", as 
for instance, in the "brusselator" model for the so called "chemical clocks" [4,5]. 

Limit cycles arise also in systems with noise. Noise or fluctuations, that are present 
everywhere, have been generally considered as a factor that destroys order. However, a 
wealth of investigations on nonlinear physics during the last decades have shown numerous 
examples, both in zero- and higher dimensional systems, of nonequilibrium systems where 
noise plays an "ordering" role. In such cases, the transfer of concepts from equilibrium ther- 
modynamics, in order to study phenomena away from equilibrium, is not always adequate 
and many times is misleading. Some examples of such nonequilibrium phenomena are: noise 
induced unimodal-bimodal transitions in some zero dimensional models (describing either 
concentrated systems or uniforms fields) [6], shifts in critical points [7], stochastic resonance 
in zero-dimensional and extended systems [8,9], noise-delayed decay of unstable states [10], 
noise-induced spatial patterns [11], noise induced phase transitions in extended systems [12], 
etc. 

Here, we discuss an extended system described by PDE's, where noise plays a key role 
in controlling and inducing a limit cycle. The model that we analyze here is the one used in 
[13,14] to study a ratchet-like transport mechanism arising through a symmetry breaking, 
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noise-induced, nonequilibrium phase transition. In a recent paper [15] a system with a 
noise induced phase transition, based on a model that is a variant of Kuramoto's model for 
coupled phase oscillators [16]; was analyzed. In addition to the phenomenon of anomalous 
hysteresis, evidence of the existence of a limit cycle for a given parameter region; is also 
shown. 

The model we analyze consists of a system of periodically coupled nonlinear phase os- 
cillators with a multiplicative white noise. Coupled oscillators have been used to model 
systems with collective dynamics exhibiting plenty of interesting properties like equilibrium 
and nonequilibrium phase transitions, coherence, synchronization, segregation and cluster- 
ing phenomena. In this particular model a ratchet-like transport mechanism arises through 
a symmetry breaking, noise-induced, nonequilibrium phase transition [13], produced by the 
simultaneous effect of coupling between the oscillators and the presence of a multiplicative 
noise. The symmetry breaking does not arise in the absence of any of these two ingredients. 
In [13] it was also shown that the current, as a function of a load force F, produces an 
anomalous (clockwise) hysteresis cycle. Recently we have reported that, changing the mul- 
tiplicative noise intensity Q and/or the coupled constant K , a transition from anomalous 
to normal (counter-clockwise) hysteresis is produced [14]. The result was obtained exploit- 
ing a mean field approximation. The transition curve in the plane (Kq,Q), separating the 
region where the hysteresis cycle is anomalous from the one where it is normal, was clearly 
determined. 

Here we focus on the time behavior. We use a method for detecting the existence of a 
limit cycle based on the evaluation of the distance between two solutions separated by a 
(fixed) time interval [17]. In this way, we not only show the existence of a limit cycle for 
F > F c (with F c a loading threshold value), but also determine its period. We also found the 
time dependence of the probability distribution function along the cycle and calculate the 
order parameter of the model vs. t, clearly showing the limit cycle. Next, we gain insight 
into its origin through the study of the large coupling limit (K — > oo). Finally, we draw 
some conclusions. 
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II. THE MODEL, MEAN FIELD AND THE METHOD USED 



For completeness we present a brief description of our model, which is similar to the 
one used in Refs. [13] and [14]. We consider a set of globally coupled stochastic differential 
equations (to be interpreted in the sense of Stratonovich) for N degrees of freedom (phases) 
Xi(t) 

f)TJ. 1 N 

This model can be visualized (at least for some parameter values) as a set of overdamped 
interacting pendulums. The second term in Eq. (1) considers the effect of thermal fluctu- 
ations: T is the temperature of the environment and the £j(t) are additive Gaussian white 
noises with 

(&(*)> = o, (&(0&(O> = M(*-*')- (2) 

The last term in Eq. (1) represents the interaction force between the oscillators. It is 
assumed to fulfill K(x — y) — —K{y — x) and to be a periodic function of x — y with period 
L = 2tt. We adopt [13,14] 

K(x)=K smx, K >0. (3) 

The potential Ui(x,t) consists in a static part V(x) and a fluctuating one. Gaussian white 
noises r)i(t), with zero mean and variance 1, are introduced in a multiplicative way (with 
intensity Q) through a function W(x). In addition; a load force F, producing an additional 
bias, is considered 

Ui(x,t) = V(x) + W(x)y/2Qrk(t) - Fx. (4) 

In addition to the interaction K(x — y), V(x) and W(x) are also assumed to be periodic 
and, furthermore, to be symmetric: V(x) = V(— x) and W(x) = W(—x). This last aspect 
indicates that there is no built-in ratchet effect. The form we choose is [13,14] 
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V(x) = W(x) = — cosx — Acos2x. (5) 

We introduce a mean-field approximation (MFA) similar to the one used in Ref. [14]. 
The interparticle interaction term in Eq. (1) can be cast in the form 
1 N 

— J2 K{Xi - Xj) = K [d(t) sin X t - Si(t) cos X,] . (6) 
i=i 

For N — > oo, we may approximate Eq. (6) in the Curie-Weiss form, replacing Ci(t) = 
N' 1 Ej cos Xj(t) and Si(t) = N' 1 J2j amXj(t) by C m = (cosXj) and S m = (sinXj), re- 
spectively. As usual, both C m and S m should be determined by self-consistency. This 
decouples the system of stochastic differential equations (SDE) in Eq. (1) which reduces to 
essentially one Markovian SDE for the single stochastic process X{t) 

X = R(X) + S(X)ri(t), (7) 

with (hereafter, the primes will indicate derivatives with respect to x) 

R{x) = -V\x) + F- K m {x) 

— - sin x(l + K C m + AA cos x) + K S m cos x + F, (8) 

(where K m (x) = K [C m sin x — S m cosx]) and 



S( x ) = V /2{T + Q[W'(x)} 2 } = ^2{T + g[sin x + 2A sin 2x} 2 }. (9) 

The Fokker-Planck equation (FPE) associated with the SDE in Eq. (7) (in Stratonovich's 
sense) is 

d t P(x, t) = d x {-[R{x) + l -S{x) S'(x)\P(x, 0) + \d xx [S 2 (x)P(x, t)} (10) 

where P(x,t) is the probability distribution function (PDF). 

In [14] we have shown that in the so called "interaction driven regime" (IDR) -where 
the hysteretic cycle is anomalous- and for each F value, in addition to the two stationary 
stable solutions with the corresponding values of current there are other three unstable ones. 
Two of them merge with the two stable, yielding a closed curve of current vs. F. Beyond a 



critical (absolute) value of the load force F, indicated by F c , those stable solutions disappear. 
This does not happen for the "noise driven regime" (NDR) -where the hysteretic cycle is 
normal-, where for each F value, one stationary stable solution exists (for small |F| even 
two stationary stable solutions and an unstable one exist). 

It is worth remarking here that the absence of a stationary stable solution, beyond the 
critical value F c in the IDR, suggest the possibility that a limit cycle exists. Already in 
[13], in a strong coupling analysis (that is considering the limit K — > oo), it was indicated 
that for very large \F\ the probability distribution function approaches a periodic long time 
behavior. 

In order to analyze the existence of a limit cycle, we exploit a novel method used in Ref. 
[17]. It is based on the measurement of the distance between different solutions of a system 
and evaluating its evolution in time. The approach applied in Ref. [17] uses a generalization 
of the known Kullback-Leibler information function [18], which is based on the nonextensive 
thermostatistics proposed by Tsallis [19]. Within such a formalism, the exponential and 
logarithmic functions are generalized according to the following definitions [17] 

exp q {x) = [1 + (1 - q)x\ 1 ' { - 1 -^ 
T i-<? _ i 

\n q {x) = - T — . (11) 

The distance can be measured between an evolved initial condition and a known stable 
stationary solution, or between two solutions at different times (separated by a time interval 
At which is fixed along the whole calculation). In this work we choose the later. In Ref. 
[17] the following definition for the distance between two solutions of a reaction-diffusion 
equation was adopted (valid for both indicated criteria) 

Pr(x,t) 



I q (Pt + Ar, P t) = ~J Pt+Ar(x, t + At) 



ln„ 



dx, (12) 



P t+Ar (x,t + Ar)_ 

where P represent a (probability-like) distribution (necessary to use the information theory 
formalism), evaluated at t and t + At, according to the criterion that we have chosen. We 
used this definition of distance, and evaluated I q (P t+ A T , P t ), using for P the PDF obtained 
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solving the FPE Eq. (10). We adopted q = 2, as it is the value for which the sensibility 
of the method seems to be a maximum [17]. The FPE was numerically solved with a 
Runge-Kuta method, using a time step 5t = 6.25 10" 7 and a space interval Sx = 0.02944. 
We have tested that variations in both steps, 5t and 5x, produce no changes in our results. 
Remembering that C m and S m should be determined self-consistently, at each time step both 
were calculated with the modified PDF. As our initial condition we adopted one stationary 
solution for F < F c calculated as in Ref. [14]. The integral in Eq. (12) was calculated 
simultaneously. Furthermore, we also obtained v m -the particle mean velocity- 



Figure 1 shows I q (P t+ Ar, Pt) (normalized to its maximum) vs. t for A — 0.15, T = 2, 
K = 10, Q = 3 and F = 1.5 (a set of parameters for which a stationary stable solution does 
not exists: see Fig 6 in Ref. [14]). We observe that I q is a periodic function of time. This 
form seems to be typical for limit cycles as shown in Ref. [17], the period corresponding to 
the distance between peaks. In Fig. 2, for the same parameter values, we depict the PDF 
at different times along the complete cycle, where the behavior resembles a wave. In Fig. 3 
we show v m and S rn vs. t. They have a time periodic behavior, not as in the case F < F c , 
where v m 0) and S m 0) are both constants in time. We have also verified that the 
transition to the limit cycle occurs just at F c (in this case F c = 1.2). 




L/2 



(13) 



which is adopted as the order parameter like in [14]. 



III. RESULTS 



A. Numerical results 



B. Asymptotic strong coupling analysis 



In order to understand the origin of the previous results and gain some insight about 
them, we have performed an asymptotic strong coupling analysis. That is, we consider 
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Kq — > oo, P — > 5(x — x m ), hence Eq. (13) transforms into 

1 

x m = i2(x m ) + -S(x m )S'(:r m ). (14) 

A simple calculation shows 

x m = — sin x rn [l + AA cosx m ] [1 — Q cosx m — AAQ(l — 2 sin 2 x m )\ + F. (15) 
This equation can be analyzed considering an effective potential U(x m ), given by 

U(x m ) = V(x m ) - QW'\x m )/2 - Fx m , (16) 
that allows us to rewrite Eq. (14) as 

*» = (17 » 

Figure 4 shows the solution of Eq. (15), x m vs. t, for both situations: just below and above 
F c . It was observed that while for F < F c , after a transient, the solution becomes stationary, 
for F > F c it becomes oscillatory. In the first case x m is constant in time but it does not 
imply v m — because, it should be calculated with S m = sin(:r m ) ^ 0, not as in the case 
with x m . Figure 5 shows the effective potential U vs. x m for the same cases, and also for 
F — 0. It is apparent that in the first case (F < F c ) the potential has only one minimum 
while for the second one, both possible minima are washed out. The latter happens just 
when the transition to the oscillating regime occurs. It is worth remarking here that, if 
Kq — * oo, the hysteresis cycle is anomalous and closed, and a critical load force establishing 
a threshold for a limit cycle transition always exists. 



IV. CONCLUSIONS 

A wealth of papers have reported on research where, by changing a control parameter, a 
transition to a limit cycle occurs [20]. However, studies on the existence of limit cycles under 
(or induced by) the influence of noise are scarce [15,21,22]. Such an aspect was analyzed 
here, where we have studied a system of periodically coupled nonlinear oscillators with 
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multiplicative white noises, yielding a ratchet-like transport mechanism through a symmetry- 
breaking, noise-induced, nonequilibrium phase transition [13,14]. The model includes a load 
force F, used as a control parameter, so that the picture of the current vs. F shows hysteretic 
behavior. 

In [14] we have found that in the IDR the cycle is anomalous, yielding a closed curve 
current vs. F when the stationary stable solutions merge with two of the three unstable 
ones. For F > F c (force value at which a stable solution merges with an unstable one) there 
are no stationary stable solutions. Here we have shown, by analyzing the time evolution of 
the distance between different solutions, that at F = F c a transition to a limit cycle occurs. 
Such a distance shows, for F > F c , a typical periodic behavior evidencing a limit cycle [17]. 
Focusing on the analysis of the time behavior, we have shown the evolution of both the 
PDF and the current, showing in both cases the time periodicity (a time evolution of the 
PDF resembling a wave). In order to understand the origin of this transition, we have made 
a "strong coupling" limit analysis. It indicates that the minima of the effective potential 
are "washed out" as F is increased and all the stationary stable solution are removed with 
them. 

As indicated in the introduction, limit cycles balance the in- and out- energy flows - 
even when those flows are not oscillatory in time- through a system's oscillatory motion 
that equalize such flows over one period. In the present case we have found a limit cycle 
in a dynamical system described by PDE's, where the energy inflow is provided by both 
the load force F and the noise terms, while energy is lost (as the system is an overdamped 
one) proportionally to the particle's velocity. A remarkable aspect is the fact that it is the 
multiplicative noise intensity the parameter controlling the bifurcation towards the limit 
cycle. 

Summarizing, for this model (that is just one example among many possible others) we 
have found a transition towards a limit cycle induced by both, a multiplicative noise and 
a global periodic coupling. However, when the noise or coupling are not present, such a 
transition does not happen. This is a new feature of those systems showing a ratchet-like 
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transport mechanism arising through a symmetry-breaking, noise-induced, nonequilibrium 
phase transition. Also, it is another example where the presence of a multiplicative noise 
contributes to build up some form of order. 
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FIGURE CAPTIONS 



Figure 1: I q (P t+ & T , P t ) (divided by its maximum) vs. time (t) for A = 0.15, T = 2, 
K — 10, Q — 3 and F — 1.5 (for this set of parameters there is no stationary stable 
solution). 

Figure 2: PDF (P) vs. x for different time (t) following the complete cycle (starting at 
t = 2.83, and evaluated each At = 0.2444). The parameters A, T, K , Q and F as in Fig. 1. 

Figure 3: v m and S m vs. time (t). The parameters A, T,K ,Q and F as in Fig. 1. Line 
thick for v m and thin for S m . 

Figure 4: Solution of Eq. (15) x m vs. t, for both situations, just below and above 
F c = 1.2. The parameters A, T, _K" and Q as in Fig. 1. We observe that while for F < F c , 
after a transient, the solution becomes stationary, for F > F c it is oscillatory. The parameters 
are K = 10 and Q = 3. The solid line indicates the case F > F c and the dashed one the 
case F < F c . 

Figure 5: U vs. x m just below and above of F c = 1.2. Also the case F = is shown. The 
parameters are A = 0.15, 7f = 10, and Q = 3. It is observed that in the first case (F < F c ) 
the potential has at least a minimum, while for the second one both possible minima are 
washed out. The solid line indicates the case just above F c (F c = 1.2), the dotted indicates 
the case F < F c and dashed one the case F — 0. 
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FIG. 2. 



15 




16 




17 




18 



